A compendium of multi-omics data illuminating host responses to lethal human virus infections.

Human infections caused by viral pathogens trigger a complex gamut of host responses that limit disease, resolve infection, generate immunity, and contribute to severe disease or death. Here, we present experimental methods and multi-omics data capture approaches representing the global host response to infection generated from 45 individual experiments involving human viruses from the Orthomyxoviridae, Filoviridae, Flaviviridae, and Coronaviridae families. Analogous experimental designs were implemented across human or mouse host model systems, longitudinal samples were collected over defined time courses, and global multi-omics data (transcriptomics, proteomics, metabolomics, and lipidomics) were acquired by microarray, RNA sequencing, or mass spectrometry analyses. For comparison, we have included transcriptomics datasets from cells treated with type I and type II human interferon. Raw multi-omics data and metadata were deposited in public repositories, and we provide a central location linking the raw data with experimental metadata and ready-to-use, quality-controlled, statistically processed multi-omics datasets not previously available in any public repository. This compendium of infection-induced host response data for reuse will be useful for those endeavouring to understand viral disease pathophysiology and network biology.

Modeling PAH Mixture Interactions in a Human In Vitro Organotypic Respiratory Model.

One of the most significant challenges in human health risk assessment is to evaluate hazards from exposure to environmental chemical mixtures. Polycyclic aromatic hydrocarbons (PAHs) are a class of ubiquitous contaminants typically found as mixtures in gaseous and particulate phases in ambient air pollution associated with petrochemicals from Superfund sites and the burning of fossil fuels. However, little is understood about how PAHs in mixtures contribute to toxicity in lung cells. To investigate mixture interactions and component additivity from environmentally relevant PAHs, two synthetic mixtures were created from PAHs identified in passive air samplers at a legacy creosote site impacted by wildfires. The primary human bronchial epithelial cells differentiated at the air-liquid interface were treated with PAH mixtures at environmentally relevant proportions and evaluated for the differential expression of transcriptional biomarkers related to xenobiotic metabolism, oxidative stress response, barrier integrity, and DNA damage response. Component additivity was evaluated across all endpoints using two independent action (IA) models with and without the scaling of components by toxic equivalence factors. Both IA models exhibited trends that were unlike the observed mixture response and generally underestimated the toxicity across dose suggesting the potential for non-additive interactions of components. Overall, this study provides an example of the usefulness of mixture toxicity assessment with the currently available methods while demonstrating the need for more complex yet interpretable mixture response evaluation methods for environmental samples.

PM2.5 Is Insufficient to Explain Personal PAH Exposure.

To understand how chemical exposure can impact health, researchers need tools that capture the complexities of personal chemical exposure. In practice, fine particulate matter (PM2.5) air quality index (AQI) data from outdoor stationary monitors and Hazard Mapping System (HMS) smoke density data from satellites are often used as proxies for personal chemical exposure, but do not capture total chemical exposure. Silicone wristbands can quantify more individualized exposure data than stationary air monitors or smoke satellites. However, it is not understood how these proxy measurements compare to chemical data measured from wristbands. In this study, participants wore daily wristbands, carried a phone that recorded locations, and answered daily questionnaires for a 7-day period in multiple seasons. We gathered publicly available daily PM2.5 AQI data and HMS data. We analyzed wristbands for 94 organic chemicals, including 53 polycyclic aromatic hydrocarbons. Wristband chemical detections and concentrations, behavioral variables (e.g., time spent indoors), and environmental conditions (e.g., PM2.5 AQI) significantly differed between seasons. Machine learning models were fit to predict personal chemical exposure using PM2.5 AQI only, HMS only, and a multivariate feature set including PM2.5 AQI, HMS, and other environmental and behavioral information. On average, the multivariate models increased predictive accuracy by approximately 70% compared to either the AQI model or the HMS model for all chemicals modeled. This study provides evidence that PM2.5 AQI data alone or HMS data alone is insufficient to explain personal chemical exposures. Our results identify additional key predictors of personal chemical exposure.

Predicting personal PAH exposure using high dimensional questionnaire and wristband data.

BACKGROUND: Polycyclic aromatic hydrocarbons (PAHs) are a class of pervasive environmental pollutants with a variety of known health effects. While significant work has been completed to estimate personal exposure to PAHs, less has been done to identify sources of these exposures. Comprehensive characterization of reported sources of personal PAH exposure is a critical step to more easily identify individuals at risk of high levels of exposure and for developing targeted interventions based on source of exposure. OBJECTIVE: In this study, we leverage data from a New York (NY)-based birth cohort to identify personal characteristics or behaviors associated with personal PAH exposure and develop models for the prediction of PAH exposure. METHODS: We quantified 61 PAHs measured using silicone wristband samplers in association with 75 questionnaire variables from 177 pregnant individuals. We evaluated univariate associations between each compound and questionnaire variable, conducted regression tree analysis for each PAH compound and completed a principal component analysis of for each participant's entire PAH exposure profile to determine the predictors of PAH levels. RESULTS: Regression tree analyses of individual compounds and exposure mixture identified income, time spent outdoors, maternal age, country of birth, transportation type, and season as the variables most frequently predictive of exposure.

Expanding the access of wearable silicone wristbands in community-engaged research through best practices in data analysis and integration.

Wearable silicone wristbands are a rapidly growing exposure assessment technology that offer researchers the ability to study previously inaccessible cohorts and have the potential to provide a more comprehensive picture of chemical exposure within diverse communities. However, there are no established best practices for analyzing the data within a study or across multiple studies, thereby limiting impact and access of these data for larger meta-analyses. We utilize data from three studies, from over 600 wristbands worn by participants in New York City and Eugene, Oregon, to present a first-of-its-kind manuscript detailing wristband data properties. We further discuss and provide concrete examples of key areas and considerations in common statistical modeling methods where best practices must be established to enable meta-analyses and integration of data from multiple studies. Finally, we detail important and challenging aspects of machine learning, meta-analysis, and data integration that researchers will face in order to extend beyond the limited scope of individual studies focused on specific populations.

Expanding the access of wearable silicone wristbands in community-engaged research through best practices in data analysis and integration.

Wearable silicone wristbands are a rapidly growing exposure assessment technology that offer researchers the ability to study previously inaccessible cohorts and have the potential to provide a more comprehensive picture of chemical exposure within diverse communities. However, there are no established best practices for analyzing the data within a study or across multiple studies, thereby limiting impact and access of these data for larger meta-analyses. We utilize data from three studies, from over 600 wristbands worn by participants in New York City and Eugene, Oregon, to present a first-of-its-kind manuscript detailing wristband data properties. We further discuss and provide concrete examples of key areas and considerations in common statistical modeling methods where best practices must be established to enable meta-analyses and integration of data from multiple studies. Finally, we detail important and challenging aspects of machine learning, meta-analysis, and data integration that researchers will face in order to extend beyond the limited scope of individual studies focused on specific populations.

3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism.

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection.

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

Multi-omics of NET formation and correlations with CNDP1, PSPB, and L-cystine levels in severe and mild COVID-19 infections.

The detailed mechanisms of COVID-19 infection pathology remain poorly understood. To improve our understanding of SARS-CoV-2 pathology, we performed a multi-omics and correlative analysis of an immunologically naïve SARS-CoV-2 clinical cohort from blood plasma of uninfected controls, mild, and severe infections. Consistent with previous observations, severe patient populations showed an elevation of pulmonary surfactant levels. Intriguingly, mild patients showed a statistically significant elevation in the carnosine dipeptidase modifying enzyme (CNDP1). Mild and severe patient populations showed a strong elevation in the metabolite L-cystine (oxidized form of the amino acid cysteine) and enzymes with roles in glutathione metabolism. Neutrophil extracellular traps (NETs) were observed in both mild and severe populations, and NET formation was higher in severe vs. mild samples. Our correlative analysis suggests a potential protective role for CNDP1 in suppressing PSPB release from the pulmonary space whereas NET formation correlates with increased PSPB levels and disease severity. In our discussion we put forward a possible model where NET formation drives pulmonary occlusions and CNDP1 promotes antioxidation, pleiotropic immune responses, and vasodilation by accelerating histamine synthesis.

Coupling Environmental Whole Mixture Toxicity Screening with Unbiased RNA-Seq Reveals Site-Specific Biological Responses in Zebrafish.

Passive sampling device (PSD) extracts paired with developmental toxicity assays in Danio Rerio (zebrafish) are excellent sensors for whole mixture toxicity associated with the bioavailable non-polar organics at environmental sites. We expand this concept by incorporating RNA-Seq in 48-h post fertilization zebrafish statically exposed to PSD extracts from two Portland Harbor Superfund Site locations: river mile 6.5W (RM 6.5W) and river mile 7W (RM 7W). RM 6.5W contained higher concentrations of polycyclic aromatic hydrocarbons (PAHs), but the diagnostic ratios of both extracts indicated similar PAH sourcing and composition. Developmental screens determined RM 6.5W to be more toxic with the most sensitive endpoint being a "wavy" notochord malformation. Differential gene expression from exposure to both extracts was largely parallel, although more pronounced for RM 6.5W. When compared to the gene expression associated with individual chemical exposures, PSD extracts produced some gene signatures parallel to PAHs but were more closely matched by oxygenated-PAHs. Additionally, differential expression, reminiscent of the wavy notochord phenotype, was not accounted for by either class of chemical, indicating the potential of other contaminants driving mixture toxicity. These techniques offer a compelling method for non-targeted hazard characterization of whole mixtures in an in vivo vertebrate system without requiring complete chemical characterization.

The Superfund Research Program Analytics Portal: linking environmental chemical exposure to biological phenotypes.

The OSU/PNNL Superfund Research Program (SRP) represents a longstanding collaboration to quantify Polycyclic Aromatic Hydrocarbons (PAHs) at various superfund sites in the Pacific Northwest and assess their potential impact on human health. To link the chemical measurements to biological activity, we describe the use of the zebrafish as a high-throughput developmental toxicity model that provides quantitative measurements of the exposure to chemicals. Toward this end, we have linked over 150 PAHs found at Superfund sites to the effect of these same chemicals in zebrafish, creating a rich dataset that links environmental exposure to biological response. To quantify this response, we have implemented a dose-response modelling pipeline to calculate benchmark dose parameters which enable potency comparison across over 500 chemicals and 12 of the phenotypes measured in zebrafish. We provide a rich dataset for download and analysis as well as a web portal that provides public access to this dataset via an interactive web site designed to support exploration and re-use of these data by the scientific community at http://srp.pnnl.gov .

Missing data in multi-omics integration: Recent advances through artificial intelligence.

Biological systems function through complex interactions between various 'omics (biomolecules), and a more complete understanding of these systems is only possible through an integrated, multi-omic perspective. This has presented the need for the development of integration approaches that are able to capture the complex, often non-linear, interactions that define these biological systems and are adapted to the challenges of combining the heterogenous data across 'omic views. A principal challenge to multi-omic integration is missing data because all biomolecules are not measured in all samples. Due to either cost, instrument sensitivity, or other experimental factors, data for a biological sample may be missing for one or more 'omic techologies. Recent methodological developments in artificial intelligence and statistical learning have greatly facilitated the analyses of multi-omics data, however many of these techniques assume access to completely observed data. A subset of these methods incorporate mechanisms for handling partially observed samples, and these methods are the focus of this review. We describe recently developed approaches, noting their primary use cases and highlighting each method's approach to handling missing data. We additionally provide an overview of the more traditional missing data workflows and their limitations; and we discuss potential avenues for further developments as well as how the missing data issue and its current solutions may generalize beyond the multi-omics context.

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 h post-fertilization after exposure to eight retene concentrations. We aimed to identify the most sensitive transcriptomic responses and to determine whether this approach could uncover gene sets uniquely differentially expressed at concentrations which induce a phenotype. We identified a concentration-response relationship for differential gene expression in both number of differentially expressed genes (DEGs) and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-ß (TGF-ß) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module seven genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module seven genes are largely responsible for retene toxicity.

Impact of acute temperature and air pollution exposures on adult lung function: A panel study of asthmatics.

BACKGROUND: Individuals with respiratory conditions, such as asthma, are particularly susceptible to adverse health effects associated with higher levels of ambient air pollution and temperature. This study evaluates whether hourly levels of fine particulate matter (PM2.5) and dry bulb globe temperature (DBGT) are associated with the lung function of adult participants with asthma. METHODS AND FINDINGS: Global positioning system (GPS) location, respiratory function (measured as forced expiratory volume at 1 second (FEV1)), and self-reports of asthma medication usage and symptoms were collected as part of the Exposure, Location, and Lung Function (ELF) study. Hourly ambient PM2.5 and DBGT exposures were estimated by integrating air quality and temperature public records with time-activity patterns using GPS coordinates for each participant (n = 35). The relationships between acute PM2.5, DBGT, rescue bronchodilator use, and lung function collected in one week periods and over two seasons (summer/winter) were analyzed by multivariate regression, using different exposure time frames. In separate models, increasing levels in PM2.5, but not DBGT, were associated with rescue bronchodilator use. Conversely DBGT, but not PM2.5, had a significant association with FEV1. When DBGT and PM2.5 exposures were placed in the same model, the strongest association between cumulative PM2.5 exposures and the use of rescue bronchodilator was identified at the 0-24 hours (OR = 1.030; 95% CI = 1.012-1.049; p-value = 0.001) and 0-48 hours (OR = 1.030; 95% CI = 1.013-1.057; p-value = 0.001) prior to lung function measure. Conversely, DBGT exposure at 0 hours (ß = 3.257; SE = 0.879; p-value>0.001) and 0-6 hours (ß = 2.885; SE = 0.903; p-value = 0.001) hours before a reading were associated with FEV1. No significant interactions between DBGT and PM2.5 were observed for rescue bronchodilator use or FEV1. CONCLUSIONS: Short-term increases in PM2.5 were associated with increased rescue bronchodilator use, while DBGT was associated with higher lung function (i.e. FEV1). Further studies are needed to continue to elucidate the mechanisms of acute exposure to PM2.5 and DBGT on lung function in asthmatics.

The Ahr2-Dependent wfikkn1 Gene Influences Zebrafish Transcriptome, Proteome, and Behavior.

The aryl hydrocarbon receptor (AHR) is required for vertebrate development and is also activated by exogenous chemicals, including polycyclic aromatic hydrocarbons (PAHs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). AHR activation is well-understood, but roles of downstream molecular signaling events are largely unknown. From previous transcriptomics in 48 h postfertilization (hpf) zebrafish exposed to several PAHs and TCDD, we found wfikkn1 was highly coexpressed with cyp1a (marker for AHR activation). Thus, we hypothesized wfikkn1's role in AHR signaling, and showed that wfikkn1 expression was Ahr2 (zebrafish ortholog of human AHR)-dependent in developing zebrafish exposed to TCDD. To functionally characterize wfikkn1, we made a CRISPR-Cas9 mutant line with a 16-bp deletion in wfikkn1's exon, and exposed wildtype and mutants to dimethyl sulfoxide or TCDD. 48-hpf mRNA sequencing revealed over 700 genes that were differentially expressed (p < .05, log2FC > 1) between each pair of treatment combinations, suggesting an important role for wfikkn1 in altering both the 48-hpf transcriptome and TCDD-induced expression changes. Mass spectrometry-based proteomics of 48-hpf wildtype and mutants revealed 325 significant differentially expressed proteins. Functional enrichment demonstrated wfikkn1 was involved in skeletal muscle development and played a role in neurological pathways after TCDD exposure. Mutant zebrafish appeared morphologically normal but had significant behavior deficiencies at all life stages, and absence of Wfikkn1 did not significantly alter TCDD-induced behavior effects at all life stages. In conclusion, wfikkn1 did not appear to be significantly involved in TCDD's overt toxicity but is likely a necessary functional member of the AHR signaling cascade.

Evaluating predictive relationships between wristbands and urine for assessment of personal PAH exposure.

During events like the COVID-19 pandemic or a disaster, researchers may need to switch from collecting biological samples to personal exposure samplers that are easy and safe to transport and wear, such as silicone wristbands. Previous studies have demonstrated significant correlations between urine biomarker concentrations and chemical levels in wristbands. We build upon those studies and use a novel combination of descriptive statistics and supervised statistical learning to evaluate the relationship between polycyclic aromatic hydrocarbon (PAH) concentrations in silicone wristbands and hydroxy-PAH (OH-PAH) concentrations in urine. In New York City, 109 participants in a longitudinal birth cohort wore one wristband for 48 h and provided a spot urine sample at the end of the 48-hour period during their third trimester of pregnancy. We compared four PAHs with the corresponding seven OH-PAHs using descriptive statistics, a linear regression model, and a linear discriminant analysis model. Five of the seven PAH and OH-PAH pairs had significant correlations (Pearson's r = 0.35-0.64, p ≤ 0.003) and significant chi-square tests of independence for exposure categories (p ≤ 0.009). For these five comparisons, the observed PAH or OH-PAH concentration could predict the other concentration within a factor of 1.47 for 50-80% of the measurements (depending on the pair). Prediction accuracies for high exposure categories were at least 1.5 times higher compared to accuracies based on random chance. These results demonstrate that wristbands and urine provide similar PAH exposure assessment information, which is critical for environmental health researchers looking for the flexibility to switch between biological sample and wristband collection.

Atomic Force Microscopy and Infrared Nanospectroscopy of COVID-19 Spike Protein for the Quantification of Adhesion to Common Surfaces.

The COVID-19 pandemic has claimed millions of lives worldwide, sickened many more, and has resulted in severe socioeconomic consequences. As society returns to normal, understanding the spread and persistence of SARS CoV-2 on commonplace surfaces can help to mitigate future outbreaks of coronaviruses and other pathogens. We hypothesize that such an understanding can be aided by studying the binding and interaction of viral proteins with nonbiological surfaces. Here, we propose a methodology for investigating the adhesion of the SARS CoV-2 spike glycoprotein on common inorganic surfaces such as aluminum, copper, iron, silica, and ceria oxides as well as metallic gold. Quantitative adhesion was obtained from the analysis of measured forces at the nanoscale using an atomic force microscope operated under ambient conditions. Without imposing further constraints on the measurement conditions, our preliminary findings suggest that spike glycoproteins interact with similar adhesion forces across the majority of the metal oxides tested with the exception to gold, for which attraction forces ∼10 times stronger than all other materials studied were observed. Ferritin, which was used as a reference protein, was found to exhibit similar adhesion forces as SARS CoV-2 spike protein. This study results show that glycoprotein adhesion forces for similar ambient humidity, tip shape, and contact surface are nonspecific to the properties of metal oxide surfaces, which are expected to be covered by a thin water film. The findings suggest that under ambient conditions, glycoprotein adhesion to metal oxides is primarily controlled by the water capillary forces, and they depend on the surface tension of the liquid water. We discuss further strategies warranted to decipher the intricate nanoscale forces for improved quantification of the adhesion.

Gene co-expression network analysis in zebrafish reveals chemical class specific modules.

BACKGROUND: Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response. RESULTS: Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure. CONCLUSIONS: Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

Unfolded Protein Response Inhibition Reduces Middle East Respiratory Syndrome Coronavirus-Induced Acute Lung Injury.

Tissue- and cell-specific expression patterns are highly variable within and across individuals, leading to altered host responses after acute virus infection. Unraveling key tissue-specific response patterns provides novel opportunities for defining fundamental mechanisms of virus-host interaction in disease and the identification of critical tissue-specific networks for disease intervention in the lung. Currently, there are no approved therapeutics for Middle East respiratory syndrome coronavirus (MERS-CoV) patients, and little is understood about how lung cell types contribute to disease outcomes. MERS-CoV replicates equivalently in primary human lung microvascular endothelial cells (MVE) and fibroblasts (FB) and to equivalent peak titers but with slower replication kinetics in human airway epithelial cell cultures (HAE). However, only infected MVE demonstrate observable virus-induced cytopathic effect. To explore mechanisms leading to reduced MVE viability, donor-matched human lung MVE, HAE, and FB were infected, and their transcriptomes, proteomes, and lipidomes were monitored over time. Validated functional enrichment analysis demonstrated that MERS-CoV-infected MVE were dying via an unfolded protein response (UPR)-mediated apoptosis. Pharmacologic manipulation of the UPR in MERS-CoV-infected primary lung cells reduced viral titers and in male mice improved respiratory function with accompanying reductions in weight loss, pathological signatures of acute lung injury, and times to recovery. Systems biology analysis and validation studies of global kinetic transcript, protein, and lipid data sets confirmed that inhibition of host stress pathways that are differentially regulated following MERS-CoV infection of different tissue types can alleviate symptom progression to end-stage lung disease commonly seen following emerging coronavirus outbreaks. IMPORTANCE Middle East respiratory syndrome coronavirus (MERS-CoV) causes severe atypical pneumonia in infected individuals, but the underlying mechanisms of pathogenesis remain unknown. While much has been learned from the few reported autopsy cases, an in-depth understanding of the cells targeted by MERS-CoV in the human lung and their relative contribution to disease outcomes is needed. The host response in MERS-CoV-infected primary human lung microvascular endothelial (MVE) cells and fibroblasts (FB) was evaluated over time by analyzing total RNA, proteins, and lipids to determine the cellular pathways modulated postinfection. Findings revealed that MERS-CoV-infected MVE cells die via apoptotic mechanisms downstream of the unfolded protein response (UPR). Interruption of enzymatic processes within the UPR in MERS-CoV-infected male mice reduced disease symptoms, virus-induced lung injury, and time to recovery. These data suggest that the UPR plays an important role in MERS-CoV infection and may represent a host target for therapeutic intervention.

Hypergraph models of biological networks to identify genes critical to pathogenic viral response.

BACKGROUND: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems. Hypergraphs are generalizations of graphs that naturally model multi-way relationships and have shown promise in modeling systems such as protein complexes and metabolic reactions. In this paper we seek to understand how hypergraphs can more faithfully identify, and potentially predict, important genes based on complex relationships inferred from genomic expression data sets. RESULTS: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges represent significantly perturbed genes, and vertices represent individual biological samples with specific experimental conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. CONCLUSIONS: Our results demonstrate the utility of using hypergraphs to represent complex biological systems and highlight central important responses in common to a variety of highly pathogenic viruses.